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Abstract 

We study a quenched SU(2) lattice gauge theory in which, in an attempt to distinguish 
between timehke and spacehke gauge fields, the gauge ensemble {f/^} is generated 
from a 3 dimensional gauge-Higgs model, the timelike link variables being "recon- 
structed" from the Higgs fields. The resulting ensemble is used to study quenched 
quark propagation with non-zero chemical potential fi; in particular, the quark den- 
sity, chiral and superfiuid condensates, meson, baryon and gauge-fixed quark propa- 
gators are all studied as functions of fi. While it proves possible to alter the strength 
of the inter-quark interaction by changing the parameters of the dimensionally re- 
duced model, there is no evidence for any region of parameter space where quarks 
exhibit deconfined behaviour or thermodynamic observables scale as if there were a 
Fermi surface. 
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1 Introduction 



Lattice QCD at non-zero quark chemical potential ought in principle to be much 
more straightforward than the corresponding problem with non-zero temperature T. 
The reason is that fi can be introduced via a local term in the Lagrangian, whereas 
T > is imposed via a lattice of finite extent L^a = in the temporal direction. 
Simulations with varying /i can therefore be performed at fixed lattice spacing a, 
so that the renormalisation factors needed for bulk thermodynamic observables such 
as the Karsch coefficients required to extract the physical energy density need be 
calculated only once, rather than for each value of T studied. Moreover, the regime 
fi/T ^ 1 relevant for the physics of the superdense matter found in neutron star cores 
implies lattices with large Lt, so that temporal correlators can be sampled efficiently 
without recourse to anisotropic lattices requiring careful callibration. This means that 
excitations above the ground state of the system, or more generally the nature of the 
spectral density function, can be studied with relative ease. Indeed, this programme 
has been successfully carried out in certain theories such as the NJL and related 
models, where both measurement of the superfluid gap in the quasiparticle spectrum 
at ~ |1] and the identification of a phonon excitation with uj (x k in the spin-1 
meson channel [2] have proved possible. 

In practice, of course, the reason this happy state of affairs has not been exploited 
in QCD is the notorious Sign Problem associated with most Euclidean field theories 
having a non-zero density of a conserved charge. In brief, for lattice QCD with Nj 
quark fiavors described by a Lagrangian density of the form qMq, the functional 
measure det^-^M(/i) = det^-^ M*(—/i) implying that for 7^ the action is complex, 
rendering Monte Carlo importance sampling impracticable in the thermodynamic 
limit. Can one then at least perform lattice QCD simulations in the quenched Nj 
limit, i.e. study the propagation of valence quarks with fi ^ through a non- 
perturbative gluon background? 

In principle the information extracted from such an approach could be at best 
qualitative, since (unlike the case of T > 0) the gauge field ensemble {U^} can 
only respond to /i 7^ via virtual quark loops, so that in an orthodox quenched 
simulation the gluon background is that of the vacuum with zero baryon charge 
density. Nonetheless, such information might be valuable, for instance, in furnishing 
a non-perturbative definition of the Fermi surface, whose existence is assumed in most 
phenomenological treatments of dense matter. 
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In the NJL studies referred to above, the Fermi surface appears as a minimum of 
the dispersion relation uj{k) in the neighbourhood of the Fermi momentum k = kp'- for 
weakly-interacting massless quarks we expect kp ^ jj,. Even free quark propagation 
at /i 7^ reveals the presence of a Fermi surface and so is not entirely trivial. For 
QCD however, the notion of a distinguished quark momentum is not gauge invariant, 
so that the identification of regions of phase space occupied by low-energy excitations 
may be altogether more subtle in a gauge invariant formalism. 

However, as plausible as this sounds, the quenched approach has conceptual dif- 
ficulties at yU 7^ 0. In the context of a random matrix theory, Stephanov [3] showed 
that the quenched theory should be thought of as the Nf — > limit of a QCD-like 
theory with not just Nf flavors of quark g G 3 of the SU(3) gauge group but also 
with Nf flavors of conjugate quark q'^ G 3. As has been known for many years |4j, 
the extra particle content results in gauge invariant qq'^ bound states in the spectrum 
which in a strongly interacting theory can result in baryons degenerate with light 
mesons. At /i = these extra states are usually regarded as extra pions due to the 
fermion species doubling resulting from the use of the manifestly real positive measure 
detMM''' required by practical fermion algorithms. The introduction of /i 7^ distin- 
guishes baryons from mesons. A simple argument, which assumes that binding energy 
is a relatively small component of the energy density of bulk nuclear matter, predicts 
that for /i/T ^ 1 there should be an onset transition from the vacuum to a ground 
state with quark number density > once /i equals the mass of the lightest baryon 
divided by the number of quark constituents. For QCD this scale is mjv/3, where 
TTiN is the nucleon mass; for a theory with conjugate quarks the onset scale is m.,t/2. 
Only cancellations among configurations due to a complex-valued measure detM can 
ensure that the fake signal for > vanishes in the range <yU<mAr/3[5]. A 

recent analytic demonstration has been given, once again in the context of a random 
matrix model, in [6]. 

It is clear from these considerations that quenched calculations at /i 7^ can 
only be useful if the distribution of gluon fields is modified in some way to reflect 
high baryon density. For instance, a minimum requirement is that the timelike links 
f/o are no longer drawn from the same sampling distribution as spacelike links Ui, 
reflecting the breakdown of Lorentz invariance due to the preferred rest frame of the 
background quark distribution. If the gluons were modified in some way so that color 
confinement no longer holds, then the role of qq'^ excitations may not be so important 
in determining the ground state in the quark sector, and it is at least conceivable that 
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valence quark propagation in such a background may qualitatively resemble that of the 
deconfined regime of the phase diagram at /i/T ^ 1 corresponding to quark matter. 
The focus of such a study would thus be "cold quarks in hot glue" . While there are 
several lines one could take, the approach we shall explore in this study is to start with 
the 3d configurations characteristic of the deconfined phase found at T > Tc, /i/T ^ 1 
produced by the approach to hot gauge theory known as Dimensional Reduction (to 
be reviewed below in Sees. 12.11 and 12.21) . In this approach all non-static modes of the 
gauge theory (ie. those with non-zero Matsubara frequency) are integrated out leaving 
a 3d gauge-Higgs model describing the non-perturbative behaviour of the remaining 
static modes. For sufficiently large T/Tc the model coefficients are perturbatively 
calculable functions of T, /i and Nf, and the resulting effective theory can be used 
to make quantitative predictions in the quark gluon plasma phase. Our goals are 
less ambitious and more speculative; we will use the Higgs fields of the 3d theory to 
"reconstruct" the timelike gauge fields Uq and hence generate {3+l)d configurations 
of static gauge fields suitable for the study of valence quark propagation with /i 7^ 0. 
At this stage we are simply trying to generate a non-confining gluon background, and 
in no sense claim to be developing a high density effective theory. 

With gauge group SU(3) and /i 7^ 0, the dimensional reduction (DR) machinery 
yields a cubic higgs self-interaction with imaginary coefficient proportional to Nf [7j, 
which is how the 3d effective theory inherits the sign problem from {3+l)d QCD. 
Whilst there is no reason to suppose the sampling of the model's configuration space 
would be any less problematic than that of the full theory, at least the complex phase 
of each configuration is now expressed in terms of a local term in the action, making 
its calculation cheap. Moreover, since the coefficient of the complex term can be 
made arbitrarily small ad hoc, the effect of gradually introducing the sign problem, 
and the interplay of the complex phase with physical observables, could be explored. 
In this paper, however, we focus on the technically simpler case of gauge group SU(2). 
In this case the action is real for even Nf (and the cubic DR higgs self-interaction 
vanishes), permitting orthodox lattice simulations with standard algorithms, e.g. [8]. 

QCD with gauge group SU(2), often referred to as Two Color QCD or QC2D, 
has been studied at ^ ^ with lattice simulations by several groups using a variety 
of formulations and algorithms O El El [101 III Il2] • Since q and q fall in equivalent 
representations of the gauge group, so that in effect q and q'^ are identical, hadron 
multiplets contain both qq mesons and qq, qq baryons, which are degenerate at /i = 0. 
In the chiral limit the lightest hadrons are Goldstone bosons, and can be analysed 
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using chiral perturbation theory (xPT) [l3]. At leading order for fi ^ T a second 
order onset transition from vacuum to matter consisting of tightly bound diquark 
scalar bosons is predicted at exactly /io = 771^^/2. At the same point the chiral 
condensate (qq) starts to fall below its vacuum value, and a non-vanishing diquark 
condensate (qq) develops, such that (gg)^ + (qq)'^ remains constant. The diquark 
condensate spontaneously breaks U(l) baryon number symmetry, so the resulting 
ground state is superfluid. In the limit /i fio+ the matter in the ground state 
becomes arbitrarily dilute, weakly-interacting, and non-relativistic, and is a textbook 
example of Bose-Einstein condensation. This scenario was subsequently confirmed by 
simulations with staggered lattice fermions [5],|9]. More recent simulations have found 
evidence for a second transition at larger /i to a deconfined phase, as evidenced by a 
non- vanishing Polyakov loop [11] and by a fall in the topological susceptibility |12j . 
In this regime thermodynamic quantities scale according to the expectations of free 
field theory (also referred to as "Stefan-Boltzmann" (SB) scaling), namely Uq oc /i^, 
and energy density e oc fj,^ 

In this paper we build on our existing experience by exploring yU 7^ in quenched 
QC2D. In Sec. Elbelow we specify our procedure for generating {3+l)d quenched SU(2) 
gauge configurations starting from a 3d gauge-Higgs model, and review standard 
quark observables once /i 7^ 0. Our results follow in Sec. [3l First we explore the 
quark density Uq and the chiral (qq) and superfluid {qq) condensates to see how the 
equation of state responds to attempts to render the {3+l)d theory "non-confining", 
and whether the behaviour predicted by xPT can be supplanted by the SB scaling 
expected of weakly-interacting degenerate quarks. Next we turn to spectroscopy, 
calculating both "normal" and "anomalous" components of the quark propagator, 
taking advantage of the enormous gain in statistical accuracy offered by the quenched 
approach. First we study the bound state spectrum in both meson and diquark 
sectors, finding evidence for significant mixing between the sectors at large enough fi. 
Next, for the first time in a gauge theory, we present results for the quark propagator 
at ;U 7^ obtained following gauge-fixing. A range of spatial momenta are sampled in 
an attempt to map out the quasiquark dispersion relation. While we find significant 
qualitative differences between strongly and weakly-interacting quarks, no signal for 
a Fermi surface has emerged. Our conclusions follow in Sec. HI 
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2 Formulation 



2.1 Reconstructing the Fourth Dimension 

The quenched action we start from is the 3d SU(2) gauge - adjoint Higgs model given 
by Eqn. (4) of Ref.[Il]: 

S3d = +2Y,^t{VxV,) 

x,i>j ^ ^ x 

- 2kJ2 ^<'{'xUx,^'{>x+^Ul^ + A ^(2tr ((^,(^,) - 1)2, (1) 



where = \^a'Ta is a traceless hermitian 2x2 matrix representing the adjoint Higgs 
field. As is usual in a gauge-Higgs model, increasing k at sufficiently large 13 takes one 
from a "confinement" phase with small (tryx^) to a "Higgs" phase with large (inpip). 

The action ([1]) is derived by dimensional reduction from a Ad pure gauge SU(2) 
theory with non-zero temperature T . In continuum notation the original Lagrangian 
density £4 ~ [d^Ayf in terms of a gluon field with canonical mass dimension [A\ = 1. 
The dimensionally reduced Lagrangian is obtained by integrating over Euclidean time 
and discarding all non-static modes, ie. setting do = 0: 
j_ 

£3 = jj dxoC, ~ ^{D.A^y + ^{D.Aof = {D.B^f + (A$)', (2) 

with the new fields defined by Bi = T-^Ai, $ = T-^Aq, [B] = [$] = !. When 
transcribing to the 3d lattice action ([1]) we use 

Ui = exp{igaAi) = exp{igaT2 Bi); if = \ (3) 

with [U] = [(p] = 0. The 3d lattice parameters /3, k, and A are all dimensionless, with 
/3 given in terms of the Yang- Mills coupling by 

4 4 AN-r 
P=—2^^ = ^^ (4) 

where it is helpful to distinguish between the dimensionful coupling g^, of the 3d theory 
and the dimensionless coupling g of the parent hot 4d theory, related via g^ = g^T; 
in the final equality we have assumed that the 4d theory is formulated on a lattice 
with Nr time spacings. 

In the DR approach to the effective description of hot field theory, the parameters 
P, K and A are perturbatively calculable functions oig, T and quark chemical potential 
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fi, as outlined in [7]. The attractiveness of this approach is that effects of, say Nf 
flavors of massless quark or a small quark chemical potential /i can be incorporated 
in the coefficient calculation, the result always being a 3d bosonic model which is 
relatively cheap to simulate. The DR provides a good effective description when the 
non-static modes decouple; the authors of [71 [H] claim that this is valid for T ^ 2Tc, 
fi ^ 4T. The parameter choice employed in this study is discussed further below in 
SecO 

As discussed in Sec. [H the quenched approximation has long been considered to be 
inapplicable to QCD with yU 7^ 0, because it includes unphysical light baryons formed 
as bound states of g G 3 and g'^ G 3 [S, H]. The reason these qq'^ states are light and 
hence distort the physics is because at T = quenched QCD configurations are in 
a confining phase, implying spontaneous chiral symmetry breaking. The lightest gg^ 
state is degenerate with the Goldstone pion. Our goal is to study quark propagation 
through a non- confining quenched gluon background with yU 7^ 0. Since chemical 
potential couples to quarks via the timelike component of the current fiip'joip: this 
is an inherently four- dimensional problem. In order to generate such a background 
we take a 3d configuration generated by the DR simulation, motivated by the fact 
that it describes deconfining physics, and "reconstruct" the gauge field in the timelike 
direction following ([3]) via the prescription (recall uAq = aT5$ = [aKT)'iip): 

Uo = exp {^9\j^'^ = cos {(jy/Lpa'fa) + ^~^= ^^"^ idV^a^a) , (5) 

with g = ^J^. Spatial link variables Ui are taken to be time independent and identical 
to their 3a counterparts. The resulting model differs from DR in that it has a non- 
trivial electrostatic potential (ie. a spatially- varying Uq field), but more importantly 
differs from real physics with /i ^ T in that there are no non-static (ie. Oq 7^ 0) 
modeil. It is not derivable from QCD in any systematic way. 

It is legitimate to ask whether excluding non-static modes from consideration is 
too severe an approximation to be physically reasonable. It is possible to examine 
this issue for degenerate quark matter interacting weakly via gluon exchange |16j . 
The loop integral in the self-consistent equation for the color-superconducting gap 
A is dominated by small angle scattering between quarks at antipodal points of the 
Fermi sphere, diverging as In^ where 6 ~ qo/f^, go being the timelike momentum 
of the exchanged gluon. The divergence must be cut off by some physical screening 



^Even the static modes in dense matter with a sharp Fermi surface may exhibit osciUatory 
behaviour, known as Friedel oscillations J^, which are absent from the DR approach. 
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mechanism: for electric gluons this is Debye screening, resulting in a ~ 9 below 
which the interaction is effectively point-like; for magnetic gluons the mechanism is 

2 - 1 

Landau damping resulting in 9min ^ 9^Qo / ■ In the static limit go — ^ magnetic 
gluons are thus unscreened in perturbation theory, and the resulting gap equation 
yields A ~ /iexp(— 37r^/v^(yf) [T7j. The effect of neglecting non-static modes such as 
electric gluons can be shown to affect the pre-exponential factor, but not the scaling of 
the gap with the coupling strength. Neglect of non-static modes thus seems defensible 
in a weak-coupling approach to quark matter; beyond perturbation theory we have 
little to guide us. 

A conceptual point worth stressing is that due to the different ways spacelike and 
timelike links are treated, there is no simple relation between spatial and temporal 
lattice spacings. In principle one could determine the relation empirically by com- 
parison of correlators in different directions, and a^/aj callibrated by matching, say, 
to the ratio of the physical (and T, /^-dependent) Debye screening mass to a hadron 
mass measured at yU = 0. We make no attempt to follow this prescription in this 
exploratory study: all masses and energy scales are compared with a reference scale 
chosen to be the pion mass at yU = 0, calculated in dimensionless units as m.,^at. This 
prevents us for the time being from a direct comparison with condensates such as the 
physically observable quark density, expressed in dimensionless units Note, 
however, that the dimensionless ratio n/T can be expressed without ambiguity using 
lattice variables as Nt-ii, where A^,- is the temporal extent of the lattice and jl = fiat 
is the dimensionless lattice chemical potential introduced in Eqn. ffT6l) below. 

To illustrate both the benefits and potential pitfalls of this approach it is help- 
ful to consider rectangular Wilson loops. The spatial Wilson loop W{ri,r2) clearly 
inherits the behaviour of the 3d theory; it is known to decay with an area law |18j . 
Confinement in the {3+l)d theory, however, is governed by the temporal loop W{r, t) 
given by 

W{r,t) = ^trP(0,0;f,0)f/*(f)pt(o,t;f,t)(f/t(0))*. (6) 

Here P is the path-ordered product of spacelike links connecting to r; because the 
configuration is static these can be written P(0, 0;r, 0) = P(0, t;r, t) = P. Using 
( I4ll5l) we deduce 

Ac 

Now expand the exponentials as power series in g, noting that (2^9)^ = (/^aV^al = Iv^pl, 
and tr^?" = for odd n, implying that odd terms vanish. The lowest non-trivial term 
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is 0{g^): 

- ^(l^or + Wr?) + ^tr(P^oPV.)- (8) 

To proceed we make some simplifying approximations. For any 7^ 0, r > 0, the 
second term of ([8]) is proportional to lPqLp\P^^] in the strong coupling limit /5 — 0, 
strong link field fluctuations should make its expectation value very small, and hence 
we neglect all terms of this form. Secondly, when integrating over fluctuations of 
the <y9 we assume translational invariance and neglect all non-Gaussian fluctuations 
(which should be valid as A ^ 0), so that (ly^oP^lv'rP^) ~ (|</?p(P+9)). This also 
requires physical excitations of the "id theory to be massive so that {^po^Pr) can be 
neglected. It is then straightforward to show that all r-dependence drops out and the 
expectation value of the 0(^^") term is (-l)"(2^t(|(^|))2"/(2 x (2n)!). Hence 

{Wir,t))^cos'Cgtm). (9) 

As required, the temporal Wilson loop is positive, well-behaved for gt\ip\ <^ 1, and 
clearly does not obey an area law. The reconstructed {3+l)d theory is thus non- 
confining. Note that the neglected PipP^cp terms always contribute to the series 
coefficients with the opposite sign (due to an extra factor of i^), so that their inclu- 
sion would in effect make (W) decay less rapidly. Viewed as a function of Euclidean 
time (W) decays with a negative curvature; this is inconsistent with a transfer ma- 
trix with positive definite spectrum, which implies decays of the form ^2- Cie~^^^. 
Unsurprisingly, therefore, the {3+l)d model violates unitarity. 

2.2 Parameter Choice for the Dimensionally Reduced Model 

The DR procedure for SU(N) gauge theory with Nf flavours of fermions has been 
studied to two loops in [I4]. The effective theory obtained is a three dimensional 
SU(N) adjoint Higgs theory, with the scalars corresponding to the electric gauge 
potential Aq in the unreduced theory; the case with chemical potential /i 7^ is dealt 
with for example in Ref. [7j. 

The continuum action of the SU(2) adjoint Higgs model is given by jH] 

S = jdl'x l^tr (Fi.Fij) + tr (A$A$) + ^'tr ($$) + A3 [tr ($$)]' | , (10) 

where Fij = diAj - djAi + ig3[Ai,Aj] and = + ig3[Ai,^. The physical 
properties of the theory are fixed by the two dimensionless ratios 

A3 mo , , 

^ = 7l> y = -T- (11) 

ys ^3 
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In the framework of DR these parameters are completely determined by and T of 
the original four dimensional gauge theory. Since g'^ is a running coupling its value is 
fixed by the renormalisation scale Aj;jg. Choosing the renormalisation scale as in [19] 
and expressing Aj^ through the critical temperature = 1.23(ll)Aj^ as measured 
on the lattice [20], it is possible to show [I9] that for Nj = 



9l 



10.7668 



0.3636 



X 



y[x) 



97r^x Att"^ 



+ — + Oix). (12) 



ln(8.3T/Tj ' ln(6.6T/Tc; 

With these equations, specifying T/Tc completely fixes the parameters x,y of the 
reduced model. Note that the values of x and y corresponding to are about 
Xc ~ 0.19 and yc ~ 0.14. The model exhibits two phases: a symmetric phase with full 
SU(2)-type confinement and a Higgs phase with residual U(l)-type confinement |21j . 

Simulations of the three dimensional SU(2) adjoint Higgs model showed that 
this theory has a phase diagram with continuously connected Higgs and confinement 
phases, which are partially separated by a line of first order phase transitions [T9l[22] . 
The discretized form of the action ( fTOj) in terms of a rescaled lattice field (p is simply 
our action ([T]). The 3d lattice parameters f3, k and A are all dimensionless. The pa- 
rameters of the continuum and lattice theory are related up to two loops by a set of 
equations [23] . 



X 



+ 



167r2 



p = 


4 

2 ' 


-3- 


2xk\ 
P ) 


(20a;- 


- lOx^) 



S/3 

477 

In 



^ 5 
1 + -X 
4 



3P_ 



+ 0.09 +8.7 + 11.6a; 



(13) 



where S = 3.17591. Due to the theory's superrenormalisability, these perturbative 
relations are exact in the continuum limit and, based on experience, are accurate for 
all /3 > 6 [2]. In Tab. [T]we report some sample lattice parameter sets obtained by 
calculating x and y according to Eqs. ( fT2l) and fixing the ratio T/Tc, then we have 
calculated k and A using Eqs. (fT3|l . 

For gauge group SU(3) the DR theory is very similar to Eq. (fTOl) [7]: 



S = jd^x jitr^ + tr[A, Ao]^ + m^trA^ 



2\2 



(14) 



where Fij = diAj — djAi + ig[Ai, Aj], Di = di + ig-^Ai, Fij, Ai, and Aq are all traceless 
3x3 Hermitian matrices {Aq = A^Ta, etc), and g^ and A3 are the gauge and scalar 
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X 


y 




A 


1 


0.193 


0.142 


0.3637812 


0.0028379 


2 


0.141 


0.185 


0.3620027 


0.0020531 


3 


0.122 


0.210 


0.3612335 


0.0017689 


5 


0.104 


0.242 


0.3603992 


0.0015009 


10 


0.087 


0.285 


0.3594500 


0.0012490 


100 


0.056 


0.427 


0.3569257 


0.0007927 



Table 1: The lattice parameters for (3 = 9.0. 



coupling constants, with mass dimension one. The physical properties of the effective 
theory, also in this case, are determined by the dimensionless ratios 

. = ^. y='^3(El^, (15) 



where Jx^ is the MS dimensional regularization scale in 3d. 

Here we note that application of DR to QCD with Nf quark flavors and non- 
zero chemical potential /i results in a complex term ig^^^tiAQ [7j, so that the 
DR theory inherits the complex action of the parent theory (for gauge group SU(2) 
tr^Q = and the action is real). In principle a term of this form could be incorporated 
in our approach, the advantage being that a "Sign Problem" could be introduced 
with arbitrarily small coefficient, and the response of observables to small phase 
fluctuations of the measure assessed. This gradualist approach is of course impossible 
in full QCD simulations. 

In applications of the DR action ([1]) to hot gauge theory, the parameter /? is used to 
control the physical lattice spacing in units of T~^, and k, and A are tuned according to 
relations ( JT2l [T3l) to specify T. A complication [Tl:| is that for the {(3, n, A} parameter 
set appropriate for hot QCD the 3d model has its true ground state in the "Higgs" 
phase, whereas the perturbative continuation to QCD f|T2|) requires it to be in the 
"confining" phas^. Fortunately the "confining" phase appears to be metastable, so 
that simulations started off with small (tryjy?) can be used to yield physically relevant 
results. Once the model is extended to Ad using ([5]) the physical meaning of the lattice 
parameters is no longer clear. In this exploratory study we hold (3, A fixed and, with 
the exception of Figs. [9] and [TO] below, restrict our attention to two values of k. 



^In this context the terms "confining" and "Higgs" refer to 3d dynamics, corresponding to the 
behaviour of spatial Wilson loops in a {3+l)d model. They have no bearing on the behaviour cither 
of the hot gauge theory before DR or the reconstructed M model studied here. 
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2.3 Introducing Quarks with fi 

Henceforth we use M configurations {U^} generated as outlined above as input 
in quenclied studies using staggered fermions in tfie fundamental representation of 
SU(2), having action S1+S2 with 

Si = IY1 X.[e'''^°U^.,.X.+o - e-^'^"Ul,^^x.-o] + mj^ X.X. = xMx (16) 



2 

and 



52 = 5^ ^Xr^2X. + (17) 



As explained above, the lattice parameter fi is related to the physical quark chemi- 
cal potential via fl = ^at. Henceforth we will ignore the distinction between a^, at 
when quoting values for parameters such as quark mass m and diquark source j; in 
consequence no attempt will be made to compare these quantities with a physical 
scale. 

The simplest observables to discuss are the chiral condensate and quark density: 

1 / dg-\\ 1 / dg-^ 

where (. . .) denotes a quenched average using the action ([1]), and Q is the Gor'kov 
propagator defined below in fl^^ . The normalisation, somewhat arbitrary in a quenched 



iSi 



two per 
is intro- 



simulation, is chosen so that in the limit /i — 00 saturates at a value o 
lattice site. The term 5*2 comprising gauge invariant diquark source term 
duced in order to discuss the possibility of diquark condensation. The Pauli matrices 
T2 act on SU(2) color indices. While 5*1 is invariant under the global U(l) rotation 
X I— > e*"x! X ^ S2 is not. Defining source strengths j± = j ± J, we obtain 

diquark condensates 

If we choose sources such that = = j, j_ =0, then the condensate forms in the 
"+" channel, and a massless Goldstone pole develops in the "-" channel as j — as 
confirmed by the Ward Identity 

J](gg_(0)gg-(^)) = (20) 



■^Note that a gauge invariant diquark source does not exist for gauge group SU(3). Physically, 
this means diquark condensation forces superfluidity for SU(2), but superconductivity for SU(3). 
Technically, it means diquark condensation is more difficult to address by lattice methods for SU(3). 
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Figure 1: vs. j from data taken on 8^ x 128 at ma = 0.05, /2 = for /3 = 9.0, 
K = 0.362007, A = 0.0020531. 



The numerical implementation of all these observables is identical to that for the NJL 
model described in [1]. 

3 Numerical Results 

We have chosen j3 = 9.0 sufficiently close to the continuum limit for the DR formalism 
to be trustworthy, and start with a point with T = 2Tc according to f|T2|) : from Tabled] 
this corresponds to n = 0.3620027, A = 0.0020531. For this parameter set, we have 
explored system volumes 8^ x Lj, with Lf = 16, 32, 64 and 128. For the most part 
in this section we present data taken with Lf = 32. However, it is important to 
make a precise determination of the pion mass at = 0. Fig. [T] shows data for 
the pion mass m-Tr as a function of diquark source j taken at yU = 0. The data are 
extrapolated to j = using a xPT-inspired form [T3] mj^^j) = m^(0)(l + bj^)^ (the 
bare quark mass ma = 0.05, so that the pion remains massive as j ^0). We conclude 
"^TrC-tO = 0) = 0.230(3), consistent with a measurement made exactly at yU = j = 0: 
m^^at = 0.2321(1). Since this scale is not too dissimilar to L~^, we have repeated the 
measurement on 16^ x 64, where we find mj^at = 0.2368(3). The systematic error 
due to finite volume is significant, but small enough at 2% to be acceptable for this 
exploratory study. 
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3.1 Equation of State 

In Figs. [2] and [3] we plot quark density Uq and chiral condensate {qq) as functions 
of /i for various j, this time on a 8^ x 32 lattice, with < /i/T < 8. The bare 
quark mass throughout this study was set to ma = 0.05. To confirm freedom from 
finite volume effects we also performed test simulations with Lg = 4, 8, 16, 32 and 
Lf = 16,32,64. There is a transition at /i ~ 0.12, becoming more abrupt as j — >■ 0, 
from a phase with Uq = 0, (qq) constant to one in which Uq increases approximately 
linearly with /i and (qq) oc This is in complete accordance with the scenario 

described by xPT in which as /i increases at T ^ there is a transition at fic = 
from the vacuum to a weakly-interacting Bose gas formed from scalar diquarks (Cf. 
Figs. 4 and 5 of Ref. [I3]). The diquarks are supposed to Bose-condense to form a 
superfluid condensate; on a finite system this must be checked at j 7^ using the 
{11+) observable of (fT9|) . Fig. H] shows a compilation of {qq+) data as a function of j 
for /i rising from zero up to /i = 0.25, the condensate increasing monotonically with 
yU. To determine the nature of the ground state an extrapolation j — > is needed. 
We have used a cubic polynomial for data with 0.02 < ja < 0.1, which may result 
in some systematic uncertainty in the immediate neigbourhood of the transition, but 
Fig. [5] confirms that once again there is an abrupt change of behaviour in the order 
parameter at /i ~ m,r/2, and that the high-/i phase is superfluid. 

Next we explored a parameter set corresponding to a smaller scalar "stiffness" by 
changing to k = 0.1 while keeping ma = 0.05 - naively following Table [1] suggests 
this corresponds to a huge value of T/Tc, ie. taking us further into the deconfined 
phase of the hot gauge theory. Of course, whether DR-based concepts remain valid 
for the reconstructed theory must be addressed empirically. This time we used a 
volume 8^ x 64 for the bulk observables and at /i = determined the pion mass 
m^at = 0.1377(1) on 8^ x 128, and m^at = 0.1423(4) on 16^ x 64, showing that 
the finite volume error is now roughly 3%. Figs. [6] and [7] show respectively Ug and 
{qq) as functions of /i for < /i/T ^ 12 in the same format as previously. It is 
noteworthy that for /i > m^/2 nq{fi) is numerically very similar to the values found 
at K = 0.3620027, whereas for /i < m.,^/2 the chiral condensate (qq) is significantly 
smaller, indicative of a weaker quark - anti-quark binding at this smaller k. As before, 
there is a clear discontinuity in the observables' behaviour at fic — m.,^/2, and the 
general picture is qualitatively very similar, suggesting that the xPT scenario is still 
applicable. Diquark binding is now also much weaker, however, as shown by the 
data of Fig. [HI Note that the source values j are greater than those of Fig. [H but the 
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Figure 6: Ug vs. fi for various j at j3 = 9.0, k = 0.1, A = 0.0020531. 



0.20 



0.16 

<qq> 

0.12 



0.08 



El \ 



i ' ' 

lattice 8 x64 

am=0.05 

p=9.0 

K=0.1 

;i=0.0020531 



GOj=0.00 0-0Hi 
□□j=0.02 ^[ 
<X>j=0.04 



0.00 



0.05 



0.10 

1^ 



0.15 



0.20 



Figure 7: {qq) vs. /i for various j using the same {/?, k. A}. 



17 



0.10- 



<qq+> 

0.05 



0.00 



1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 






- oeoe^ 




)0eeeee00eo^ 

lattice 8^64 




GOj=0.02 
□□j=0.04 
□□j=0.06 


am=0.05 
B-9.0 

K=0.1 

^1-0.0020531 

1 . , , , 1 , , , , 


1 , , , 


1 , , , , 



0.00 



0.05 



0.10 

1^ 



0.15 



0.20 



Figure 8: {qq+) vs. fi for various j using the same k, A}. 



corresponding are smaller in magnitude. Since the curvature of the data as a 

function of j is quite pronounced, a reliable extrapolation j — > is impracticable on 
this system size. Indeed, inspection by eye would suggest a linear extrapolation would 
yield {qq+) — for all /i; however, the ja = 0.02 data do manifest some significant 
variation at yU ^ ^12^^/2, suggestive that weak symmetry-breaking persists for fi > fi^- 
It is disappointing that we have found no qualitative change in physics as the 
parameters are varied - recall that the xPT model which describes the results rea- 
sonably well is based on the assumption of confinement, or at least on the presence of 
very tightly bound diquark states in the spectrum. To explore the parameter space 
more widely we focussed on a single observable, Uq, and scanned the (k. A) plane on 
8^ X 16 at five different values of /i with (3 = 9.0, ma = 0.05 and ja = 0.01. The results 
are summarised in Figs. [9] and [TOl Fig. [9] shows that except for A = 0.1 the results for 
fixed n are practically independent of k (shown by the overlapping symbols) and of A 
(shown by the horizontal lines). Fig. [TO] shows that increases linearly with /i over 
a wide region of parameter space, as it does for /x > /Zc in Figs. [2f6l This approximate 
linear behaviour is once again a prediction of xPT [HI [13], and is to be contrasted 
with the Uq oc /i^ behaviour expected of a deconfined theory where baryons can be 
identified with degenerate quark states occupying a Fermi sphere of radius kp ^ fi- 
The absence of this scaling is a further reason to conclude that the reconstructed 
model does not describe deconfined physics. 
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3.2 Bosonic Spectrum 

Next we report on bound state spectroscopy. Conventionally, qq bound states are 
referred to as mesons, and qq, qq as diquark baryons and anti-baryons respectively. In 
a medium with spontaneously broken baryon number symmetry, however, excitations 
need not have a well-defined baryon number. An appropriate set of states to look at, 
together with gauge-invariant interpolating operators expressed in terms of staggered 
fermion fields, is as follows: 

pion Xx^xXx] (21) 
scalar XxXx] (22) 

higgs ^{xi^2Xx + XxT2tx) = QQ+x] (23) 
goldstone ^{xir2Xx - Xxr2tx) = QQ-x- (24) 

Here the phase Ex = (— l)^o+^i+^2+a;3^ g^j^^j ^]-^g Pa^jj matrix T2 acts on color indices. 
Pion and scalar states are related via the U(l)e global symmetry x ^ X ^ 

^g«ae_ Analogous to chiral symmetry for continuum spinors, this is an exact symmetry 
of the action (fT6|) in the limit m 0. In a phase with spontaneously broken chiral 
symmetry, the scalar is massive, and the pion a Goldstone mode, becoming massless 
as m — i> 0. Similarly, "higgs" and "goldstone" diquark states are related via the 
U(1)b baryon number rotation x ^ ^^^Xi X ^ X^"*^, an exact symmetry of fll6|17l) 
in the limit j — > 0. In a superfiuid phase with 7^ 0, the higgs is massive, and 

the goldstone massless in the limit j — > 0. 

The boson correlators are constructed from the Gor'kov propagator [T] 



fxy 







( ^xy 


N \ 




•^^2 / xy 


\ ^xy 


Axy J 



(25) 



where the 2 x 2 (in color space) matrices ~ {XxXy) aiid A ~ {XxXy) are known 
as the normal and anomalous parts respectively. On a finite volume A = for 
j = 0; limj^o liniy->oo ^ 7^ signals particle-hole mixing resulting from the break- 
down of U(l) B symmetry, and hence superfiuidity. Due to SU(2) symmetries the only 
independent components of Q are ReA^n = A^ and Imy4i2 = A and their barred coun- 
terparts, so that in practice all information from a single source can be extracted at 
yU 7^ using just two matrix inversions per configuration. Each of the correlators C{t) 
constructed from the forms ( !2Tll24l) receives contributions from diagrams containing 
both A^ and A-tjpe propagators; by construction, however, they remain symmetric 
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Figure 11: Timeslice propagators of various bosonic bound states at /i = 0.25 on a 
8^ X 128 lattice at /? = 9.0, k = 0.3620027, A = 0.0020531, ma = 0.05, ja = 0.02. 

under t —t even once fi ^ 0. Note also that in this quenched treatment contribu- 
tions to the mesons from disconnected loops and to the diquarks from disconnected 
A loops are neglected. 

Fig. [11] shows all four timeslice correlators evaluated on a 8^ x 128 lattice with 
parameter set (3 = 9.0, k = 0.3620027 and A = 0.0020531. The chemical potential 
/i = 0.25, ie. above the critical /ic required to enter the superfiuid phase. The quark 
bare mass ma = 0.05 and diquark source ja = ja = 0.02. All four channels yield 
clear signals for single particle bound states, the higgs being the noisiest. 

Like any meson constructed from staggered fermions, the correlators in principle 
describe two states and must be fitted using the form 

C{t) = A\e-^' + e-'"(^*-*)] + fi[e-^* + (-l)*e-^(^'-*)]. (26) 

where m and M denote the masses of states with opposite parities. Fits to ( !26ll . 
where positive, are shown by solid lines in Fig. [TTl In most cases we find M ^ m; 
however for /i > /i^ the pion correlator has a distinct "saw-tooth" shape, and in fact 
the fit yields > M^i, where vr denotes the usual pseudoscalar pion, and 61 a state 
of opposite parity, which must therefore be scalar (note that the logarithmic scale 
requires Fig. [TT]to plot |C7r(t)|). 

In Fig. [12] we plot m^r and M^i against fi, and in Fig.[T3]the corresponding spectrum 
for all the states (12111241) . Fig. [T2] also shows m.„ and Mhi as measured on a 16'^ x 64 



21 



0.25 



0.20 



0.15 



0.10- 



O 



@ 



V 



O 



O 



lattice 8x8x8x128 

j=0.02 

am=0.05 

p=9.0 

K=0.3620027 
X^O.0020531 



O 



V 



O O 
V V 



o 



V ^ V 



V 



0.0 



0.1 



0.2 0.3 



0.4 



0.5 
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system at some representative values of /i. Even for the lightest measured mass, 
M^^{jl = 0.9) ~ 10 < Lg] while volume effects are statistically significant, they 
have no impact on the qualitative trends we now discuss. First note that all states 
are approximately degenerate at = 0. The equality of pion, higgs and goldstone 
correlators is guaranteed by SU(2) symmetry at /i = [8J, but the degeneracy of the 
scalar in the chirally-broken vacuum suggested by Fig. [3] can only result of 

meson-diquark mixing due to j ^ 0. Next, note that the pion mass remains constant 
for /i < fic, where as reviewed above it is a pseudo- Goldstone boson associated with 
chiral symmetry breaking, and then falls once the superfluid phase is entered. This is 
in accordance with the predictions of xPT for the so-called "P^" state of a theory with 
Dyson index /3o = 4 [13], and has also been observed in simulations with dynamical 
staggered quarks in the fundamental representation of SU(2) [9]. Most of the other 
states, including the bl in Fig. [121 show a much steeper decrease with /i for /i < /ic, 
followed by a gentle rise to a plateau at ma ^ 0.13 in the superfluid phase /i > fic, 
precisely that expected of the goldstone state expected in the superfluid phase with 
diquark source j ^ (Cf. the "Q/" state shown in Fig. 3 of [13]). The exception 
is the higgs, which rises more steeply to become the heaviest state at large fi. We 
conclude (a) the breaking of degeneracy between higgs (l23l) and goldstone states 
is clear supplementary evidence for the breaking of U(1)b symmetry in the superfluid 
phase; (b) all states with = 0+ including the bl but except the higgs have some 
projection onto the Goldstone state, regardless of whether the original interpolating 
operator is mesonic or baryonic. It would be interesting to study this phenomenon 
as j is varied. 

3.3 Fermionic Spectrum 

We have also studied the fermion spectrum, often in the context of condensed mat- 
ter called the quasiparticle spectrum. Since the Gor'kov propagator Q is not gauge 
invariant we have to specify a gauge fixing procedure. A feature of the quenched ap- 
proach is that it permits large statistics to be accumulated with relatively little CPU 
effort. This has enabled us for the first time in a gauge theory context to study Q at 
/i 7^ 0, by helping to overcome the sampling problems associated with gauge fixing. 
We have experimented with two gauge choices: Unitary gauge {p ip' = (0,0,{p'^), 
which is implemented before the reconstruction of the 4th dimension, and is unique 
up to a Z2 factor, specified by demanding ip'^ > 0; and Coulomb gauge, implemented 
by maximising tr([4_j + f/^_j •), in an attempt to make the gauge fields as smooth 
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Coulomb gauge 


xVdof 


Unitary gauge 


xVdof 


En 


0.1244(7) 


202 


0.1651(17) 


182 


Ea 


0.0674(13) 


3.4 


0.0696(15) 


5.9 



Table 2: Fitted mass values to the data of Figs. [Ml and [TSl for = 0. 

as possible and hence improve the signal-to-noise ratio. 

In Figs. [13] and [15] we plot respectively the normal and anomalous fermion times- 
lice propagators on a 32x8^x64 lattice at (3 = 9.0, k = 0.1, A = 0.0020531, ma = 0.05, 
ja = 0.0.2 and ft = 0.3, the last value chosen to ensure > fic- These plots result 
from an analysis of 800 independent configurations. As discussed in Ref. [1], the 
properties of staggered lattice fermions are such that in a phase with approximate 
chiral symmetry, as expected for fi > fi^ from Fig. [T] the numerically important con- 
tributions to N{t) are from t odd, and to A{t) from t even, and only these points are 
plotted. Most of the data was obtained using Coulomb gauge and differing values of 
the momentum k^a G {0, j^, . . . , |}, though the k^ = data taken in unitary gauge 
are shown for comparison. Two features to note are that the Coulomb data is roughly 
twice as large as the unitary data reflecting an enhanced signal as anticipated above, 
and that there is little variation with k^. 

In the NJL model the quasiparticle propagator can be successfully fitted using the 
forms [1] 

N{t) = Pe"^^* + Qe-^^^^'-^\ (27) 
A{t) = i?[e-^^* - e-^^^^'~% (28) 

where for 7^ there is no reason to expect P = Q, but for a well-defined quasiparticle 
state the equality Ejy = Ea should hold. Note that the anomalous amplitude i? 7^ is 
evidence for baryon number symmetry breaking, often in a condensed matter context 
called "particle-hole mixing" . Fits to the data of Figs. [TJ] and [TJ] are shown in 
Table [21 and it should be noted that only the anomalous channel fits produced an 
acceptable x^- Some insight can be gained from Fig. [161 which compares the Coulomb 
gauge normal and anomalous propagators at A;^ = on a logarithmic scale; while the 
fit ( l28j) for A{t) looks plausible, the normal channel never settles to a well-defined 
quasiparticle pole, even with Lt = 64. Moreover only the anomalous channel shows 
any evidence of gauge independence in Tab. [2l The value of Ea obtained is very close 
to m7r/2, indicating that at this value of k the pion is a weakly bound state. Another 
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Figure 14: Normal component of the quark propagator N{t) for various momenta k 
for fi > /ic- 
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Figure 15: Anomalous component of the quark propagator A{t) for various momenta 
for /i > fie- 
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Figure 16: Comparison of the = data of Figs. [H] and [I5] on a logarithmic scale. 



striking feature of the data is the approximate forwards-backwards symmetry of N{t), 
implying P Q. 

Fig. [I7| shows the dispersion relations EA{kx) for data taken on a 32 x 8^ x 64 
lattice at the parameter values shown (we were unable to obtain satisfactory fits to 
( [27|) to extract Ei<^ for ^ 0). It confirms that the quasiparticle excitation energies 
are A;- independent. This should be contrasted with the findings of [1], where a lattice 
study of the NJL model using identical formalism found E{k) exhibiting a pronounced 
minimum at ~ fc^ ~ yU, the Fermi momentum, evidence that the NJL model has 
a well-defined Fermi surface with kp ^ fi. One motivation for our study was to 
investigate to what extent the concept of a Fermi surface, which is not strictly gauge 
invariant, can be put on a firm empirical footing in a gauge theory. Fig. [T7l shows no 
evidence for a Fermi surface. 

The nature of the quasiparticle excitation is clarified a little at the other parameter 
set studied, namely k = 0.3620027 (with all other parameters unchanged). In this 
case our results show no evidence for any well-defined spin-| state in either normal 
or anomalous channels; as a result of confinement the excitation spectrum of the 
model seems to be saturated by the tightly bound spin-0 states of Fig. [131 Fig. [TSl 
shows a close up of A(t) for various /i values, showing the presence of an oscillatory 
component whose amplitude initially grows with fi (the jl = 0.3 points overlay those 
from fl = 0.2), but whose wavelength is roughly //-independent. The origin of the 
oscillation could possibly be associated with the non-unitarity of the model discussed 
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Figure 17: Dispersion relations Ejq{k) and EA{k). 
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Figure 19: E and T vs. jj for k = 0.3620027 in both normal and anomalous channels. 
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Figure 20: y/WTT^ vs. n for K = 0.3620027 in both normal and anomalous channels. 



in Sec. 12.11 but is most likely a manifestation of independent spin-| excitations being 
ill-defined due to confinement. We have fitted the k = 0.3620027 data to the forms 

N{t) = Pe-^^' cos(r^t + 0) + Qe-^^(^'-*) cos(rjv(Li - t) + 0), (29) 
A{t) = R[e-^^' cos{TAt + 0) - e'^^^^'-'^ cos{TA{Lt -t) + 0)], (30) 

where we interpret E as the energy and F as the width of a quasiparticle excitation. 
The results are plotted in Fig. [191 Most of the results are taken in Coulomb gauge, but 
at /i = 0.3 data from unitary gauge is also available. The results' most striking feature 
is their independence of /i, with F of the same order of magnitude as E. An interesting 
systematic effect is that E^ > Ea while F^r < F^, which has motivated us in Fig.[20]to 
plot \fW^^rV^ vs. 11. The figure has the same vertical scale as the previous two, and 
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it is clear that the disparity between normal and anomalous channels is significantly 
reduced. Inspection of the p, — 0.3 data also shows that the gauge dependence of this 
result is 0(20%) at worst. Numerically, y/E^~+T^ > m-,^, indicating strong quark - 
anti-quark binding, due to the persistence of confinement at this value of k. Therefore 
we can interpret the effect of confinement as rotating the quasiparticle pole into the 
complex plane, the rotation angle being larger in the anomalous channel than in the 
normal one. It would be interesting if this feature could be reproduced by analytic 
methods such as self-consistent solution of Schwinger-Dyson equations. 

4 Conclusion 

Our attempt to alter the nature of the gluon background by changing the parameters 
of the 3d DR gauge-Higgs model has been a partial success, in that in going from 
K = 0.3620027 to K = 0.1 the strength of the binding between quarks weakens signif- 
icantly. The main evidence for this claim comes from the spectrum; at the smaller k 
studied quasiquarks appear to be well-defined independent degrees of freedom, whose 
excitation energy is to fair precision half that of the gauge-invariant pion state. At 
K, = 0.3620027 by contrast, the quasiquark propagator exhibits a pole at complex k, 
and the resulting estimates for both energy and width of the excitation exceed the 
pion mass. However, in neither case is there evidence for significant departure of 
Ug, {qq) and {qq) from the behaviour predicted by xPT, so that even if quarks are 
important degrees of freedom at k = 0.1, there is no evidence for the formation of a 
degenerate system signalled by SB scahng. Moreover, our attempts to measure the 
quasiquark dispersion relation E{k) have been unsuccessful; while there is evidence 
for the gauge- invariance of the minimum energy or "gap" at k = 0.1, no /c-dependence 
is observed. This is, of course, entirely consistent with the well-known fact that the 
only gauge invariant feature of a propagator is the location of its poles. We conclude 
that identification of a Fermi surface, supposing one existed, presents a technical 
challenge in any gauge theory simulation. 

It is clear that our simplistic treatment, despite being manifestly gauge invariant 
and inclusive of non-pointlike interactions, has missed some essential component of 
the physics of high density. If we had first attempted the quenching programme using 
gauge group SU(3), the departures from the theoretical expectation of SB scaling in 
the large-// limit could have been ascribed cither to unexpected non-perturbative 
effects, or to some subtle cancellation due to the Sign Problem entirely missing from 
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the quenched approach. For SU(2), however, orthodox simulations of full QC2D, in 
which the Sign Problem is absent, yield evidence for deconfinement and SB scaling at 
large /i [HI [12] , which our approach misses completely. We can of course speculate on 
which important features of the gluon background at high quark density are absent; 
possible candidates are non-static modes, and modes with \k\ ^ 2kF ^ 2/i [15]. Sadly 
though, it appears to remain the case that despite its "unreasonable effectiveness" in 
virtually all other aspects of lattice QCD, the quenched approximation has nothing 
useful to tell us about the physics of high quark density. 
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